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Abstract 

A stochastic model for pure-jump diffusion (the compound renewal process) can be used as a 
zero-order approximation and as a phenomenological description of tick-by-tick price fluctuations. 
This leads to an exact and explicit general formula for the martingale price of a European call 
option. A complete derivation of this result is presented by means of elementary probabilistic 
tools. 
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I. INTRODUCTION 



Market microstructure and high-frequency trading are fields ex per iencing an increase of 
interest within financial institutions and academic scholars (see Therefore, models 

for tick-by-tick financial fiuctuations, that mi ght have been considered just a curiosity more 



than a decade ago (see the seminal pap er by 10||), are becoming more and more important 



for practical purposes (see 26|] and 17| for early studies of the model presented below). The 



literature on high-frequency price modelling up to the middle of the last decade is reviewed 



in Il2l and in 



15|- 



Along with interest in high-frequency trading comes the need for high-frequency hedging. 
In this work, we address the problem of determining the price for an intra-day European 
option written on a share traded in a stock exchange. It is assumed that the derivative 
position is opened at a time t after the start of continuous trading with maturity at a 
time Tm before the end of continuous trading on the very same day. The fiuctuations of 
the share price S{t) can be modeled by a suitable cadlag (i.e. right-continuous with left 
limit) pure-jump process. One might be tempted to use a compound Poisson process for 
the logarithm of the price X{t) = \og{S (t) / S (0)) and immediately apply the result of jl^ 
with the coefficients of drift and diffusion set to zero. For a vanishing risk-free interest rate 
(which is a reasonable assumption for intra-day data, see the discussion below), this would 
lead to the following formula for the plain- vanilla option price C (t) 

^e-^(TM-t)^ (MTM -t))" ^^ ^ ^^^^^ (1) 
n=0 

where A is the activity of the Poisson process for trades, K is the strike price, /i and are, 
respectively, the expected value and the variance of the log-price jumps which are assumed 
to be normally distributed. One further has that 

a(5(0), K, fi, a^) = iV(t/i,„)S(0) - N{d2,n)K, (2) 

where 

N{u) = —= / dve-""'/^ (3) 
v27r J^oo 

is the standard normal cumulative distribution function and, finally 

log{S{0)/K)+n{f, + ay2) 

di,n = 1= , (4) 

Vna 
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2,n 



(5) 



Whereas Merton's jump diffusion model is still the object of active research and is dis- 
cussed in a recent research paper by it has several unrealistic features. One of them 
is tha.t uncoiiditional inter-trade durations do not follow the exponential distribution (see 

Vlarkov models (described in ISj) can take this fact into ac- 
25|. In the following, we shall consider intra-day European 




count as shown by [27| and 



options written on semi-Markov pure jump models which are compound renewal processes 
25 1 and Related papers are j20| and |4|. 



20|, the focus is on option prices 



fsee also 



In 



for derivatives written on compound Poisson processes and in the presence on non-vanishing 
risk-free interest rate, whereas extends Levy option prices to the semi-Markov case by 
developing suitable approximations. Finally, a recent paper by [28| considers Laplace trans- 
form methods to deal with order and trade renewal flows in an agent-based model where 
the trade counting process is not necessarily Poisson. 

The model proposed here has several distinctive advantages; for instance, each model 
entity has an immediate and clear translation into a microstructural quantity (there are no 
hidden variables of any kind); moreover, it is very parsimonious (leading to a many sound 
mathematical results); finally, the model is very flexible. 

In what follows, section |Tl] will present the model, section IIIII will be devoted to pricing 
European options, and section IIVI will report a discussion of the results. We would like to 
stress that the derivations and computations below will involve only elementary probabilistic 
methods. 



II. AN ELEMENTARY BUT COMPREHENSIVE MODEL FOR TICK-BY-TICK 
PRICE FLUCTUATIONS 

Let S{t) denote the price of an asset at time t and let X(t) = log(S'(t)/5'(0)) be the 
corresponding logarithmic price (or log-price), where S{0) = S{t = 0) will be assumed to 
be S{0) = 1 without loss of generality. We can take as 5(0) the opening price of the asset 
after the opening auction and before the beginning of continuous trading in a stock market. 
After continuous trading begins, trades will take place at specific epochs {Tj}^£^ where M 
denotes the total number of trades within the day. We shall further assume that Tq = 0. The 
trading epochs can be seen as a point process (see ^]). Our first assumption is that these 
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epochs are a renewal process (see [6|]), meaning that the inter-trade durations Ji =Ti — Ti^i 
are independent and identically distributed (i.i.d.) positive random variables. Notice that, 
for the n-th epoch, one has 

n 

T. = 5^J„ (6) 

even if the i.i.d. hypothesis is not satisfied. If the durations {Ji}'^^ are i.i.d. random vari- 
ables, there is a simple and convenient relationship between the distribution of durations and 
the distribution of epochs. Indeed, given two independent random variables U and V with 
respective cumulative distribution functions Fjjiu) = F{U < u) and Fy{v) = F(y < v), 
one can compute the cumulative distribution function F\y{w) of their sum W = U + V, 
which turns out to be the the measure convolution (a.k.a. Lebesgue-Stieltjes convolution) 
of the two distributions for U and V (see [2]). Throughout this paper, we shall use the 
so-callled indicator-function method, a procedure allowing to write probabilities as expec- 
tations. Therefore, it is instructive to show in detail how the distribution of the sum of 
two independent random variables can be derived, even if this is a well-know result. The 
first step is to notice that the joint cumulative distribution function Fuy{u,v) is given by 
Fu,v{u,v) = Fu{u)Fy{y) as a consequence of independence. The second step is to recall 
that the probability of an event A is given by the expected value of the indicator function 
Iai namely W{A) = E(Jyi) and that the indicator function of the intersection of two events 
A and B is the product of the parent indicator functions, i.e. we have Iahb = IaIb- The 
last step is to notice that the event {W < w} is equivalent to {U G M} fl {V < w — U}. 
Therefore, one has the following chain of equalities 

Fw{w) =nW<w) = ¥. {I{w<u,}) = E {I{umhy<u,-u}) 

= [ [ dFu,viu,v)= I I dFu{u)dFy{v) 

J uGR J v<w—u J u£M. J v<w~'U 

= [ dFu{u) f dFy{v)= f Fv{w-u)dFu{u). (7) 

JweM Jv<w-u JueR 

To denote the convolution, which is an operation symmetric in U and V, we can introduce 
the symbol -k: 



Fw{w) = / Fv{w ~ u)dFu{u) = / Fu{w — v)dFv{v) 
JueR JveR 

= Fu* Fv{w) = Fv* Fu{w). (8) 
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This formula holds true also if U and V are positive random variables. In such a case one 
has that Fuiu) = for u < and that Fyiy) = for t> < (notice that -Fc/(0) and -FV(O) 
may be positive). Then equation ^ becomes 

Fw{w)= Fv{w - u)dFu{u) = Fu{w -v)dFv{v). (9) 
Jo Jo 

Now, let Fj{x) denote the cumulative distribution function of the duration, i.e. Fj{t) = 

P(J < t); moreover, let -Ft„(^) denote the cumulative distribution function of the n-th epoch, 

i.e. FT„{t) = P(T„ < t). Then, FT„(t) is given by the ra-fold convolution of Fj{t), that is by 

FT„{t) = Fr{t). (10) 

Equation fllU]) can be proved by means of the iterated application of equation (Q. 

A price S(Ti) corresponds to each trading epoch Tj. Let Yi = \og{S (Ti) / S (Ti^i)) represent 
the tick-by-tick logarithmic return, then the log-price X{t) is given by 

Nit) 

Xit) = Y,Y,, (11) 

i=l 

where the counting process N{t) is defined as 

N{t) = max{n : T„ < t}, (12) 

and counts the number of trades since the beginning of continuous trading. The relationship 
between the log-price and the price is 

Nit) 



5(t) = e^W=e^-^» = ne''^ (13) 



,JV(t) 

1=1 

,Nit) 



We shall further assume that {Fjl^jj^ is a sequence of i.i.d. random variables such that 
Y = K{Yi) < oo. Let J-'t denote the natural filtration of the process S{t) up to time t, this 
being the a- field generated by the random variables Ti, . . . , T/v(t) and Fi, . . . , Y^it)- 

In general, with the above hypotheses, S{t) given by equation f lT3|) is not a martingale. 
In fact one has for s < t 

/Nit) \ Nis) / Nit) 

E(5(t)i J-,) = E n ^ri-^. = n ^""^ n ^^-^^ 

yi=l J i=l \i=Nis)+l 

Nit) 

= S{s) n E (e^O ' (14) 

i=Nis)+l 



and the martingale condition is statisfied only if for every s, t such that s < t, one has 

N{t) 

n E(e^O = l; (15) 

i=N{s)+l 

this is the case if E(e^') = 1. However, one can always find an equivalent martingale measure 
(e.m.m.). One can replace Yi in equation f ITT]) with Yi — a defining the following processes, 
a modified log-price process 

N(t) 

X{t) = Y,{Y,-a), (16) 

1=1 

as well as the corresponding modified price process 

^(t)=e^W. (17) 
Now, if a = log(E(e^')), one has that S{t) is a martingale. In fact, one can write 

nS{t)\:F,) = S{s) W E(e^-'^) =5(5). (18) 

i=N{s)+l 

Armed with this e.m.m., it is possible to move on and price options written on the process 
defined above using the martingale method. 

III. MARTINGALE OPTION PRICING 

For an intra-day time horizon, we can safely assume that the risk-free interest rate is 
zero; even if such a return rate were ry = 10% on a yearly time horizon, meaning that the 
institution or government issuing this instrument is close to default (so that, it would not 
be so riskless, after all) or that the inflation rate is quite high, the interest rate for one day 
would be ~ 1/(10 ' 200) = 5 • 10^'* (200 is the typical number of working days in a year) 
and this number has still to be divided by 8 (number of trading hours) and then by 3600, if 
the goal is approximating the rate at the time scale of one second. This gives 1.7- 10~®. 
On the other hand, typical tick-by-tick returns in a stock exchange are larger than the tick 
divided by the price of the share. Even if we assume that the share is worth 100 monetary 
units, with a 1/100 tick size (the minimum price difference allowed), we shall have a return 
r larger than 1 ■ 10~^ and much larger than r^; therefore, it is safe to assume a vanishing 
risk-free interest rate. 
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We shall focus on the price of an intra-day European call option assuming that the 
position is taken at a time t coinciding with or close to the beginning of the day and that it 
is closed at a later fixed time (the maturity) within the same day, which we shall denote by 
Tm (not to be confused with the epochs Tj, in general Tm is not an epoch). Notice that the 
condition t < Tm must always be fulfilled. 

Let C{S{Tm)) represent the pay-off oi a European call option at maturity. For instance, 
given the strike price K, the pay-off of a plain- vanilla European option is C{S{Tm)) = 
max(0, S(Tm) —K). Then, the option price C(t) at a time t < Tm is given by the discounted 
conditional expected value of the pay-off at maturity with respect to the e.m.m., that is 

C{t)=e^^'-^-%{C{S{TM))\J't), (19) 
where r is the risk- free interest rate. In our case r = 0, so that equation (fT9|) simplifies to 

C{t)=E~^{C{S{TM))\J't)- (20) 
In order to evaluate equation fl20l) . we consider two cases: 

1. t coincides with a renewal epoch; 

2. t does not coincide with a renewal epoch, 

with the second case being the only realistic one, but the first one is discussed in the recent 
literature as a starting point for developing approximations as in \^\. 

If the option price is evaluated from a renewal epoch, we can assume that t = without 
loss of generality and the option price is given by the following integral 

^ POO ^ 

CiO) = Eg(C(5(TM))|J-o) = / C(«)dFg(^^^)(n), (21) 

Jo 

where F^^j.^^-^{u) is the cumulative distribution function of the random variable S(Tm)- In 
order to obtain this quantity, we can first define 

^„ = ne^-i°s(^(^"')); (22) 

i=l 

Sn is the product of i.i.d. random variables and its cumulative distribution function is 
the n-fold Mellin convolution of Fyiu), the common cumulative distribution function of 



16|. We shall write 



Yi = Yi — log(E(e^*)); the Mellin transform is discussed in 29|] and in 

^5„H=^r"H- (23) 
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Since the number of trades from to Tm can be an arbitrary integer, by purely probabilistic 
arguments, one can show that F^^j,^^^{u) is given by 

oo 

^§(Tm)(«) = Y.^iN(^M) = (24) 

n=0 

as a consequence of the mutual independence of tick-by-tick log-returns and inter-trade 
durations. Notice that the zero-fold Mellin convolution is a cumulative distribution function 
which is for u = and 1 for u > 0. To see that this is the case, consider equation ( 12T|1 
when it is known that n = 0. Then S(Tm) = S{0) = 1 and the payoff is C{S{0)) = C(l). 
However, the probability ¥{N{Tm) = 0) of the event A^(Tj\/) = decreases with increasing 
Tm and the contribution to the conditional expectation f l2T|) is P{N{Tm) = O)C(l). In order 
to use equation ( |24|) . we still need to compute the probabilities of the events {N(Tm) = n}. 
This can be again done by means of the indicator-function method. In fact, one has that 

{N{Tm) =n} = {T„ < Tm} n {T„+i > Tm}. (25) 

Therefore, the following chain of equalities holds true 

F{N{Tm) = n) = n{Tn < Tm} n {T„+i > Tm}) = E (/{t„<t„}/{t„+,>t„}) 



= E (/{T„<T„}/{J„+.>T„-T„}) = / / dFriu)dFjiw) 

Jo JTm-u 

/•Tm 

= / {l-Fj{TM-u))dFr{u). (26) 
Jo 

In the general case in which t is a generic observation time not coinciding with a renewal 
epoch, things become trickier, even if we are using a simplified and stylized model. At time 
t, both the price S{t) and the number of trades N{t) = nt are known. We can consider 
the random variable AX(t,TAf) = X{Tm) — X{t) = \og{S (Tm) / S (t)) . If S(t) is used as 
numeraire (that is if we set S(t) = 1), Equation ( 12T1) modifies to 

/oo ^ 

where the cumulative distribution function F~^ (u) is given by 

oo 
n=0 

Again, as in the case of equation (!2T1) . this equation can be justified by purely probabilistic 
arguments. However, one has to compute the conditional probability P(A^(TAf) — N{t) = 
n\N{t) = rtt). As derived in 13], this is given by 
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n^iTu) - N{t) = n\N{t) =nt)= I F^NiTu) -N(t + u)=n- l)dFj,^„^{u), (29) 







where F{N(TM) — N(t+u) = n — 1) is given by equation fl26|) with Tm replaced by TM — (t+u) 
and Fj^ (u) = P(j7j,„j < u) is the cumulative distribution function of the residual life-time 
at time t conditioned on the fact that there were Ut trades up to time t which we denote by 
J^t.rif The residual life time is the time interval from t to the next renewal epoch TN{t)+i- 
As discussed in [l^, its distribution crucially depends on what is known of the previous 
history. In our specific case, as anticipated above, it is meaningful to assume that we do 
know the total number of trades up to time t, as this is usually public information. Before 
deriving Fj^^^^{u), it is important to discuss the meaning of equation fl29l) . The right-hand 
side contains the probability of having n — 1 trades between the renewal epoch t + u and 
maturity Tm- Since the value u of the residual life time J7i,„t is not known, this probability 
must be convolved with the probability of the event {J^t,nt = u}. It turns out that even the 
cumulative distribution function Fj^^^^ (u) can be found by direct elementary probabilistic 
tools without using Laplace-tranform methods. We can see that the event {J^t,nt ^ can 
be described in term of a conditional event (see 8|) 

{Jt,n, <u} = {T„,+i - t < u\N{t) = rit}. (30) 

Equation (l30l) can be written in terms of epochs using ( |25l) 

{Jt,m <u} = {T„,+i -t<u\{Tr,,<t}n {T„,+i >t}}. (31) 

One can now use the definition of conditional probability and the indicator-function method 
to compute Fj^ ,^ (u) directly. First of all, one can write 



Fj,,Ju) = P(J-i,„, <u) = P(T„,+i - t < ^z|{T„, <t}n {T„,+i > t}) 

F{{Tn,+l -t<u}r\ {Tn, <t}n {Tn, + l > t}) 



(32) 



P({T„, < t} n {T„,+i > t}) 
and the denominator is already given by equation (126|1 . meaning that one has 

P({T„, <t}n {T„,,+i > t}) = [\l - Fjit - w))dFT'{w). (33) 
In order to compute the numerator, one can use the following equality between events 
{Tn,+i-t < u}r\{Tr,, < t}n{T,,+i >t} = {T,, < t}n{t-Tn, < J„,+i < t + ^-T^J, (34) 



and obtain that 



P({T„,+i -t<u}n {T„, <t}n {T„,+i > t}) = 

n{Tn, <t}n{t- Tn, < Jn,+1 < t + U - Tn,}) = 

ft PU+t — W 

^{hTr.t<t}ht-Tnt<Jn, + l<t+U-T„,}) = dFx^^ (w) dFj (v) 

Jo Jt~w 



t PU+t — W 

/ dF*/^iw)dFjiv) 

Jt-w 



[\fj{u + t-w)- Fjit - w))dFT^{w). (35) 
Jo 

Combining equations fl33l) and fl35l) . one finally gets from equation fl32|) 



Equation fl5B]) is the last ingredient needed to determine the option price in the general case 
(I27|) . Finally, note that equation (127|) yields equation ([T]) when J ~ exp(A) and F ~ -^(z^, cr^) 
(see yj], chapter 7). 

IV. DISCUSSION AND OUTLOOK 

Equation (1271) is our main results and gives a martigale price for intraday European 
options when assuming the parsimonious model of section [Tll we have been able to explic- 
itly derive all the terms in that equation by repeated application of the indicator-function 
method. Even if such an equation may seem cumbersome, we already showed that it can 
be used in practice in [l^. A more detailed numerical analysis will be the subject of future 
research. 

However, some assumptions in section [Til are unrealistic even if they ensure analytical 
tractability. For example, it is assumed that the durations { Jj}^]^ and the tick-by-tick log- 
returns {Yj^fLy are i.i.d. random variables and that they are mutually independent. In |9|]. 



3, 



221 and 



181], as well as in many other empirical papers on financial econometrics (see 
3|), it is shown that this is not the case. There is heteroscedasticity and there is dependence 
between the activity and the volatility. Suitable mixture models based on the compound 



Poisson processes can take all that into account as discussed in 24], but models using 
heteroscedastic procesess subordinated to Hawkes processes could be a viable alternative as 
well (see 2l| and references therein). 
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This paper is the crowning achievement of an activity on modelhng ultra-high frequency 
financial data by means of continuous time random walks that started back in 1998. Con- 
tinuous time random walks is the name that physicists use for compound renewal processes, 
even if some authors reserve this name to more general processes with finite or infinite mem- 
ory subordinated to counting processes. As briefiy discussed above, these processes allow 
the derivation of many non-trivial analytical results, but they are not general enough to 
take into account all the features of high-frequency financial data. In 1998, the idea was to 
use these processes for intra-day option pricing, but only in 2011, with the results published 
in [l^, it became possible to present martingale option pricing in the simple way based on 
renewal theory outlined in this paper. 
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